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nomenclatures 


fA(4>)J = periodic state matrix 
c d 888 external damping 
Cf ** mac hined clearance of journal bearing 
C ‘J' (iJ ~ x ’ y ) = damping of journal bearing 


c tj' 0. J“x,y) = dimensionless dam ping ( = 


external damping ratio f j 
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UteJL. INTtVnOMAXLV 


*(,) or m - periodic function, 0<g(t)<l and it has the form of Eq.5 for small crack. 

Jk c « complex stiffness of cracked shaft 

k t <■ stiffness of the uncrackcd shaft 

* , k - stiffness of cracked shaft in { and r, directions 

A*/ Ak, - the largest stiffness change in { and n directions caused by crack 
AK ( - stiffness change ratio in { (the crack) direction ^ - j, ^ ^ 

.... ( Ak \ 

A K n - stiffness change ratio in n (the cross) direction ^ ^ J 


AX, 


AX, 


Ak t +Ak, 


Afe, - Afe, 


^ (lj= Xi y) = stiffness of journal bearing 


c Jfc 


IF / 


jjr (i,/ 3 * xjO 35 dimensionless stiffness ^ 
mt , w , - mass lumped at bearing station and rotor mid-span 

m, - preload factor (- 1 -^-, where c» is the assembled clearance) 

r i j «... f _ OflZQk ( A V | where « is the viscosity of lubricant 

S = Sommerfeld number ^ ncre#i 

D is the bearing diameter, X is the bearing radius and L is the bearing width 

. . ( s . ( _ tHL ' D . L ( A^ \ it is a design parameter. 

So = fixed Sommcrfcld number ^ ^ \ w ) 2rc W \c f ) J 

t * time 
T ™ period 

( e -A 

U - unbalance parameter ^ g J 
W m bearing load 

W - gravity parameter ( - f . where 5, - it represents the elasticity of the 

f \ C p • 

shaft. 

x y « deflection in Cartesian coordinates 

X,Y m dimensionless deflection ( “ ~ > J~) 

. . #» 

2 s 2 +iy, complex deflection “ (s + iij)e 
Z = X+iY 

z= x-fy 

Z, * steady state response . 

AZ - perturbation deflection about the steady state response ZJ - AX + IY) ) 

/ _ \ 


mass ratio 
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P crack angle, it is the angle between crack and unbalance 
[T(J)]= transition matrix 

* body fixed rotating coordinates, { is in the direction of crack 
iif ■ Floquct eigenvalue 
\lkf | » modulus of/ijr 
t ■■ at, dimensionless time 
<p = initial unbalance angle 
<D = at + <p + p 
a = rotating speed 

= rigid pin-pin critical speed 
Q = speed ratio ( ■» — 

\ (U 

r 



(-fcn 


(' 


t 

f 


. d d 
) — -j- or — 

OT <M> 

] = square matrix 
} “ column matrix 


INTRODUCTION 

Fatigue cracking of rotor shaft may cause a catastrophic damage to rotating machineries, an example was 
given by Jack and Patterson (1976). So that a detailed investigation into the behavior of a cracked rotor shaft 
is very important for diagnosing and preventing rotor cracks. 

There are two major research field in cracked rotor analysis, one is the modeling and another is the de- 
lecting of the rotor crack. Gasch (1976) considered in body-fixed rotating coordinate the stiffness change due 
to open and closed cracks and solved the equations of motion of a cracked Jcffcott rotor by an analogue com- 
puter. Mays et al (1984) gave an approximate method of estimating the reduction of a section diameter re- 
quired to model a crack. Nelson and Nataraj (1986) analysed the dynamics of a rotor system with a cracked 
shaft by the finite element method; Collins et al (1991) considered a rotating Timoshenko shaft with a single 
transvers crack and analysed the free vibrations and the responses to a single axial impulse and periodic axial 
impulses. Gasch (1988) et al analysed the dynamic behaviour of a JefTcott rotor with a cracked hollow shaft 
and Bcmasconi (1986) considered the vibrations of a stepped shaft by using singularity functions. Muszynska 
(1982) developed a solution for a gaping crack and concluded that the increase in the primary vibration is 
more than that in the subharmonic vibrations. Papadopoulos and Dimarogonas (1980) pointed out that the 
informations of subharmonic resonances and the frequency shifting are important for crack identification. 
Imam et al (1989) presented a very successful on-line crack diagnosis method which can detect cracks of the 
order or 1 to 2 percent of the shaft diameter deep. Most of the early research results on cracked rotor have 
been summarized in the book by Dimarogonas and Paipctis (1983) and more literatures on the dynamics of 
cracked rotors can be found in a reviewing paper by Waucr (1990). 

Although more than 100 papers have been published in this field, there are still many problems to be 
solved. Up till now, only few papers, for example Tamura (1988) and Gasch (1992), analysed the stability of a 
cracked rotor supported on rigid bearings and it was found that the unstable zones caused by rotor crack are 


% 
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m a n A qj 

near the speed ratio — -7. =:. 7 ,- and the unstable zone spreading out from — - 2 is the 

CO ^ l * J “ c 

broadest, where w e is the first rigid criUcal speed. Up till now, no papers about the stability of a cracked flexi- 
ble rotor supported on fluid film bearings have been published, but the stability and dynamic characteristics of 
this kind of system arc very important as dangerous cracks were mainly reported in rotors running in fluid film 
bearings. 


EQUATIONS OF MOTION AND PERIODIC 
STATE MATRIX FOR STABILITY ANALYSIS 


Considering a simple flexible JerTcott rotor with a cracked shaft supported on journal bearings (see Fig.l), 

the equations of motion of the system can be written as: 

. , „ s «« + *) 

m/J , + c,z, + *® * 

m * ~ * {z » ~ Z P“ m >K ~ U » + ,/ r ) 

where m 4 and m d are the mass lumped at bearing station and rotor mid— span, i d and z* are the complex 
deflections of disk and bearing centers; a> is the rotating speed; * e [(z w - *» the stiffness coefficient which 

is a function of z r z h and changes with time t due to the openning and closing of the rotor crack and, f x and 
jr ^ fluid film forces of journal bearing. There are 

f z w ” x a + t *e 


(1) Stiffness Forces 

Supposing that the vibration deflection is small enough and the influence of deflection difference z 4 -z k on 
shaft stiffness k t can be omitted, then the shaft stiffness k c changes only with time t and 

M Z *- Z *), , 1 *M') (2) 

Supposing further that the weight deflection is dominant, then fc c (/) changes periodically with time and 
can be written in { and ij directions as: 

r fc ,-k - g(t)Ak 

J { ' (3) 

where #(/) is a periodic function of time and 0<g(t)<l. Ak ( and A k„ are constant and represent the 
stiffness change caused by rotor crack. There is 

z - x + iy - ({ + iif)e ** 

where G> «■ cot+<p+P, z is z d or z h . The stiffness force k c (t) • z can be written as: 


1 1 i* 

k'(t)' z-^[(k t + k n )’ z + (k ( -*,)•« ] 

-[it, -ig(0(Ak { + Afc f )]*z-|g(/XAk t -Afc e 2 * (4) 

For small crack, the openning and closing or the crack can be represented by a rectangular function g(t), and 
g(t) can be written in a Fourier series of <I> as: 

g(t) — /?(<!>) - ^ + cos< ^ — 5 — ••• j (5) 

where <l> = a)t+<p+p. 
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(2) Linearised Fluid Film Forces and Coefficients 

The vibration response of the system can be split to 


£ 




+ Az. 


* '*# ' 

wherex* andz* are the steady state response of the system, Az„ and Az* are the perturbaUon deflections 

For small dynamic perturbation, the nonlinear fluid film forces f x and f, can be linearised about the 
steady state response (x^, y u ) as: 


K]- &]*[-: •;][“:] * [:: 


( 6 ) 


The nondimensional forms or f and / can be written as 


The nondimensional stiffness and damping coefficients K {j and C tj (tj- x,y) are 


C C (0 

jr « __ l p r p 

tj fy K lj • C tj ” jy C IJ 


and they can be obtained from the book by Someya (1989) and the book by Vance (1988). These coefficients 
are only function, of Sommerfeld number 5 or eccentricity ratio for a given ret of bearing 

parameters, where the Sommerfeld number is 


S- 


fuoDLf Jt 
InW \ 


s-y 


( 8 ) 

As the rotating speed w and the bearing parameters are mixed up in the Sommerfeld number S it is diffi- 
cult to separate in the stability analysis the influences of rotating speed and bearing parameters on the 
nondimensional stiffness and damping coefficients. In order to solve this problem, the fixed Sommerfeld num- 
ber S e is introduced and it is defined as 

• \ <o ) 2nfV \ Cf ) (9) 

That is, the rotating speed or in Sommerfeld number S (see Eq.8) is replaced with the rigid pin-pin critical 

speed »(-(-) So is thus the Sommerfeld number in the rigid pin-pin critical speed and does not 

change with the rotating speed. Now the stiffness and damping coefficients are functions of the fixed 
Sommerfeld number S 0 and the rotating speed ratio 

In this paper, three tapes of journal bearing, are used for stability analysis, the 2 Anal Grooved Cylindri. 

“ 4 Lob ‘ B “ rin ‘ “ d “• 5 P * d rdtioj Pad Bearing Fig.2 show. the diagrams of them journal 

hmnngs and the stiiT.es. and damping eoefneient, of these three bearing, am taken from Someya', book 
( 1989 ). 

For the simple Jeffcott rotor considered here, there is 


W 




So that 
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w 




m c o) 

<* 9 


( 10 ) 


where W — — c -y 

CO) 

9 * 


is the gravity parameter. 


(3) Nondlmcnsional Equations of Motion of the System 

Substituting Eq.4 and Eq.6 into Eq.l, dividing both sides by m/r, co 2 and then substituting into Eq.7 and 
Eq.10, wc can get at last the nondimcnsional forms of Eq.l as 

z \ + + ^-{t 1 -^(OAJST.kz, -z,)- 

- I/KOAJr, • (Z, - z t ) * e ^ + Uc ' 




Z' 4 + -M[l - ,KZ t - Z,) - i*(/)A* 2 • (Z t - Z,) • e a *} + 

* 2«n 1 2 * 

+ ^7 +*„**. + *„A + C„AA", + C„A1“.) + 

+ K„AJr, + Jt„Al-, + C„Ar, + C, Al",)] -A 


d AJiT, + AAT a A*. -A* 

where ®nd AX, “ ^ « AJCj — ^ 

In the steady state case, there arc 

z i “ z *- z » = z ».> z '< - z 4 - °* z '» “ z '* “ °> 4 “ 4. and 4 “ 4. 

As Z^ and Z^ are also solutions of Eq.(ll) and 

Z, - Z„ + AZ,, Z, - Z„ + AZ, 

the perturbation equations of motion of the system can be got from Eq.(ll) as: 

A Z\ + A Z' 4 + ^ { t 1 - \ AK x g(t)] • (AZ , - A Z,)- 

- l -AK jf! (t)(AZ d -AZ t ) • a**} -0 

A Z\ + -h- {[1 AJe,g(/)l(AZ, - AZ; - \ AKjWAZ" - AZ,)* 2 * } + 

* %n t z ** 

[(*„A*» + *„Ar t + C„ AX\ + C X AY\) + 


2a£l 
IP (1 +2a)i 


2afl 


(ID 


f(J^A*/+ K Ay, + c„A*', + C/r,)] -0 

The above equations can be written in Cartesian coordinates as (#(/) is written in th form of £($)): 


rotor: 


n oir^/ 
Lo ij|_Ay', 


2D 

f 

n 


2D 

r 

n 


Ax 
A Y 
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bearing: 


AAT 3 sin2® 
A*, - AA' J cos2<J> 


[ AJ5T, + A^ J cos2<I» 
AA’ J sin24> 

|*(®) • 


ir- A N.o 

J \_AY,-AY t J 


( 12 ) 


'AK t + Ai5T J cos2<l> 


AAT„sin24> 


W (1 + 2a) 

E < 

2afl 2 


AAT 3 sin2© 
\K, K 

xx x f 

K K 

L y* w J 


AATj — AJf 2 cos2<X> 


A^ 

Ar 


* j 


+ 


)K=":] 
][“:)- 


(4) Periodic State Matrix of the System 
Taking the notes of 

AT, -A*, , A^-AT, , AT, 

Ar,-AAT', , A'.-A Y' t , AT 7 -A^ , X AY\ 

supposing that <p and 0 do not change with time, and renoting 0 = d/d4>, the first order state equations of 
motion of the system can be obtained from Eq.12 and then can be written in matrix form as 


•AAT, , X t ~AY k 


\XY - M(4>)){Ar} 

where the vector {Af} is defined as fA'}-(Af ) , X V ”-JY ) T 
and [>4(4>)] is the periodic state matrix and has a period of 2n. 


(14) 


[A(®) = 


where 


0 

0 

0 

0 1 

0 

0 

0 

0 

0 

0 

0 0 

1 

0 

0 

0 

0 

0 

0 o 

0 

1 

0 

0 

0 

0 

0 o 

0 

0 

1 

p t m 
n 2 

Q,(<&) 

n 2 

P,(®) 

n 2 

Q,(®> 2ft 

n 2 o 

0 

0 

0 

Q ,(<!>) 
n 2 

p 2 (<t>) 
n 2 

Q,(i>) 

n 1 

p,(«>) 

ir 0 

2/i 

n 

0 

0 

p,(a>) 

2on 2 

Q,(®) K 

2afl l 

u + p ,(®> 
2afl 2 

2aT2 

0 

c 

XX 

2 all 2 

c 

1 71— 

2afl 2 

Q,(®) 

2«n J 

P,W k^-q,^) 

2a£l 2 2afl 2 

0 

2aQ 

0 

C 

T* 

2an 2 

C 

ZL_ 

2aft 2 



P.W-l 

~ \ A*,**®) - AAr 2 g(d>)cos2® 



P 2 (®) = 1 

~ \ AAT,g((h) + ^AAT 2 g(<I))cos2® 




AK^Q) • sin 2$ 





(15) 


(16) 


The stability of the periodically time-variant system of Eq.14 can be analysed by the well known 
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Floquet's method. 


FLOQUET'S THEORETICAL AND NUMERICAL 
METHODS USED FOR STABILITY ANALYSIS 


The most straightfoward method for dealing with the stability of a periodic system is perhaps the 
Floquet's method (see eg., Caesari, 1970), which states that the knowledge of the state transition matrix over 
one period is sufficient for determining the stability of the equilibrium position f A'} = 0 of a periodic system 

fjfl'«M(*)W*> (i7) 

where [A ($)] ■» [A (d> + T)1 and T is the period. The transition matrix [T(7)] is defined by 

wr»-[r<m*(0)} 

When the transition matrix [T(7)] has been obtained, the Floquet's eigenvalues of the system can be calculated 
from the following eigen-fiinction 


|fT(7*)l-p//l|“0 (18) 

where n f is the Floquet eigenvalue and [/] is the unit matrix. 

The stability of the equilibrium position { Jf) = 0 of the periodic system (17) can be determined by the val- 
ue of the Floquet eigenvalue n, that is, if all |p/|<l, the fJf)*=0 position is stable; if one or more l/i, | > 1, the 
(JT) = 0 position is unstable and the {X} = 0 position is the critical case if one or more \fi)= 1. The physic mean- 
ing of the Floquet eigenvalue fi { is clear, it shows the factor by which a given initial amplitude increases or de- 
creases after one period of motion. 

Although the stability of the equlibrium position fA'} = 0 of system (17) can be determined from the tran- 
sition matrix [T(7)] of the system, unfortunately, there is no general analytical method for calculating [T(7)]. 
But there are several efficient numerical methods although they take too much computer time (Fridemann, 
1977). One of them is the so called Hsu's method and the essential aspects of this method are given in 
Fridemann's paper (1977), the following is the result of this method. 

The period T is divided into K intervals denoted by <b k (fc = 1,2,— JC) with 0 = ® 0 <4>, <••* <<&*“ T. The 
kth interval (G> _ ( , <b k ) Is denoted by t*. In the kth interval, a constant matrix [BJ is defined by 

[5,1 -J** M(&)](H> Get, (19) 

*k- I 


and the final approximation for the transition matrix [T(7)] at the end of one period can be calculated by 


1 r f ' [B 1 

[T(n 1 = n« P [iM* nim+ i-M 

i-l l-l l Jm\ J* ) 


(20) 


STABILITY AND STABILITY DEGREE OF A CRACKED FLEXIBLE 
ROTOR SUPPORTED ON JOURNAL BEARINGS 


Based on the Floquet's theorem and the Hsu's numerical method, the Floquet eigenvalues of Eq.15 are 
calculated numerically for speed ratio 0.2<0<5.0, crack stiffness ratio 0.0< A AT { <0.8, gravity parameter 
0.001 < W K < 10.0, fixed Sommcrfeld number 0.1 <5,,< 2.3 and mass ratio 0.05<a<0.45. The stability and the 
stability degree of a cracked flexible rotor supported on journal bearings arc analysed. At first, the following 
phraseologies are defined with Fig. 5 as an example. 

crack ridge: a protruding ridge of Floquet eigenvalue caused by rotor crack. 

crack ridge zone: the zone where the crack ridge locates, for example from A to B in Fig. 5. The position 
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and extent of the crack ridge zone is given by Eq.(21) and Eq.(22). 

crack nnatable zone : unstable zone caused by the crack, it is the zone where \n)> 1.0. It locates within the 
crack ridge zone. 

■tability degree : the degree of stability which is determined by the largest Floquet eigenvalue \n f |„ , the 
larger the I | a is, the smaller the stability degree. When | fij |„ < 1, the equilibrium position of the system is 
stable, when \p f \ m > 1, the equilibrium position of the system is unstable and the equilibrium position of the 
system is the critical case when \n f |„ ■* 1. 

diagram of stability degree : the 3-D diagram of Floquet eigenvalue. Its Z-axis is the largest Floquet 
eigenvalue \p f \ m , Y-axis is the speed ratio O and X-axis is the system parameter (the fixed Sommerfeld 
number S ot the gravity parameter W v the mass ratio a or the stiffness change ratio AK { ). 

1. 0— Level : The level of the Floquet eigenvalue \p) m * 1, it is the limit of stability. Above this leval, the sys- 
tem is unstable. 

(I) The Stability and Stability Degree of Uncrackcd Rotor Shaft 

It should be noticed that now the periodic system (17) is simplified to a time— invariant one and the 
Roquet's method is also simplified to the normal eigenvalue problem. 

Journal bearings are superior to rolling element bearings in high-speed operation. They have high damp- 
ing ability, high load-carrying capacity and low friction. In addition, silent operation and long life can be at- 
tained. Therefore, journal bearings of various types are widely used in high-speed rotating machineries. But 
an inappropriately designed journal bearing may not only increase the vibration amplitude, but bring also the 
oil whirl or oil whip, that is, a self-excited lateral vibration of rotating shaft caused by oil film in journal bear- 
ings. Generally, multilobe bearings are more stable than circular bearings. A tilting-pad bearing is always sta- 
ble, because it has no cross-coupling terms of oil film coefficients. 

The gravity parameter W t is a very heavy parameter in the stability analysis. Fig.3 shows the influences 
of gravity parameter on the stability degree. It is found that there is a valley of Floquet eigenvalue for gravity 
parameter of mid— value (0.05 < W t < 1.0) and small speed ratio (1. As O increases, this valley becomes more 
narrow and more shallow. The tilting pad bearing has the longest valley and the cylindrical bearing has the 
shortest one. The system will be the most stable if the parameters are selected in the bottom of the valley. In 
the left side (the small gravity parameter side) of the valley, there is a hill of Floquet eigenvalue for cylindrical 
and multilobe bearings. This hill is very important for system's stability as it determines the stability limit. Be- 
sides, the hill of the cylindrical bearing is much higher and wider than the multilobe bearing, that is why the cy- 
lindrical bearing is more unstable than the multilobc bearing. The tilting pad bearing is always stable as it has 
no such a hill. 

Although the multilobe bearing is generally more stable than the cylindrical bearing, it may be more 
unstable than the cylindrical bearing in the case of small fixed Sommerfeld number S and small gravity 
parameter W r 

Another heavy parameter in the stability analysis is the fixed Sommerfeld number S 0 . Fig.4 shows the in- 
fluences of S, on the stability degree in different gravity parameter W t . For multilobc bearing and for cylindri- 
cal bearing, when S t > 0.4 or W* >0.1, it is found that the threshold speed ratio of stability increases as W in- 
creases. 

The variation of stability degree with fixed Sommerfeld number S t depends largely on the value of gravi- 
ty parameter W x . For small W v the stability degree decreases as S. increases; while for large W v the variation 
of stability degree with S„ is very small. 
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(2) The Stability and Stability Degree of Cracked Rotor Shaft 


It is well known that the stability and stability degree of an uncracked flexible rotor supported on journal 
bearing change greatly for different types of journal bearings. What will happen if the rotor has a small crack? 
Of course, the crack will have influence on the stability of the uncracked system and will cause new unstable 
regions, but how large is the influence and where are these new unstable regions? These questions are impor- 
tant for detecting the crack and preventing the system From an accident. In order to answer these questions, 
the Floquet eigenvalues of Eq.15 are calculated. For the reason of simple, the cross stiffness change caused by 
crack is omitted in the Following analysis, that is A £,=*0 is assumed, because A K n is very small compared 
with AATj. 

Based on the large quantities of numerical results, it is found that the position of crack ridge zones and 
then the position of crack unstable zones have no relation with the type of journal bearings, that is, they do 
not depend on the Fixed Sommerfeld number S r The position of crack ridge zones depend only on the speed 
ratio and the stiffness change ratio for large gravity parameter ( W K > 1.0), and depend also on the mass ratio 
for small gravity parameter ( W g <0.05 ). 

Fig.5 to Fig. 7 show the positions and the extents of the crack ridge zones. When W n > 1.0, the mass ratio 


has no influence and the crack ridge zones locate always within the two lines of ft 


_2 

N 


and fl 


A*, 


than 


IN 


\ 

-^- L j, where N = 1,2, 3, 4,5,— . The widthes of the crack ridge zones are seperately equal to or less 
(see Fig.5). So that the crack unstable zones locate within 


| ( 1 - J A* { ) < fl < £ (N - 1,2,3,4,5,— ) (21) 

When W n <0.05, the crack ridge zones depend largely on mass ratio and they locate always within the two 

20 20 / i \ r\ . \ 1/2 

lines of O - and O - 1 - - AK J, where O^ - and N= 1,2,3,4,5 — . The widthes 

A K. 

of the crack ridge zones are seperately equal to or less than / O (see Fig.6). So that the crack unstable 

2 N € 

zones locate within 


20/1 \ 20 

^(l--A*J<«<^ (N - 1, 2,3,4, 5,—) (22) 

In this paper, the shaft stiffness is supposed to change linearly with crack depth, so that the width of the 
crack ridge zone increases also linearly with crack stiffness change ratio due to the reduction of shaft stiffness. 
Besides, the larger the speed ratio and the stiffness change ratio are, the wider the crack ridge zone is (see 
Fig.7). Therefore, the instability caused by crack in larger speed ratio is more dangerous and should be payed 
more attention to. 

It is known that the gravity parameter W % is a very heavy parameter in the stability analysis of uncracked 
rotor shaft and it is found in the above discussions that when W n <0.05, the mass ratio a has large influence on 
the crack ridge zone, while when W K > 1.0, a has almost no influence on the crack ridge zone. What will hap- 
pen if 0.05 < W n < 1.0? Fig.8 shows the influence of gravity parameter W % on the stability and stability degree 
(please compare Fig.8 with Fig.3). It is Found that in a region of Wg within 0.05 < W n < 1.0, the crack has al- 
most no influence on the stability degree no matter how large the crack is. This is because that in this region, 
the diagram of stability degree has either a valley or a hill. The cylindrical bearing and multilobe bearing have 
generally a valley at small speed ratio and a hill at large speed ratio, while the tilting pad bearing has always a 
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valley although this valley becomes smaller as the speed ratio increases. 

It is found that the position of crack ridge zones almost do not change with the fixed Sommerfeld number 
S t , but the extent of crack ridge zones decrease with the increasing of S t when the gravity parameter is small 
(^* < 0-05) and has almost no change or increase a little with the increasing of S„ when the gravity parameter 
is large (W t > 1.0) (see Fig.9). When comparing with the uncracked rotor case, it is found that the change of 
crack ridge zone with S, is not due to the crack itself, but due mainly to the influence of S, on the stability 
degree of uncracked system. 

Although the crack ridge zones almost do not change with the fixed Sommerfeld number S a , the crack 
unstable zones become a little bit larger as S , increases. This is also due to the change of the stability degree of 
uncracked system with S ,. 

Fig.10 shows the influences of crack stiiTness change ratio on the stability and stability degree. It is found 
again that the position of the crack ridge zones depend not on the types of journal bearings and, the crack 
unstable zones locate always within the crack ridge zones, but the stability outside the crack ridge zones de- 
pends mainly on the type or bearings used. Although the positions of crack ridge zones do not change with 
gravity parameter W M and mass ratio a, the crack unstable zones depend both on the stiffness change ratio 
and the system parameters. When W t is small, there is an extra crack unstable zones in high speed ratio due to 
the influence of mass ratio. When fV t “0.1, almost no crack ridge zones exist, because the system is now in the 
valley of the Floquet eigenvalue. So that, as it is mentioned before, the crack has almost no influence on sys- 
tem's stablity in some combinations of system parameters. As W t increases, the crack ridge zones appear 
again, but the crack ridge zone in high speed ratio is disappeared because now the mass ratio has no more in- 
fluence on the crack ridge zones. 

It is found that for tilting pad bearing, the crack ridge in small gravity parameter case is much higher and 
therefore is more unstable than in large gravity parameter case. 

CONCLUSIONS 

This paper investigates the stability and the stability degree of a flexible cracked rotor supported on dif- 
ferent kinds of journal bearings. The influences of the crack stiffness ratio A K t , the fixed Sommerfeld number 
S of the gravity parameter W t and the mass ratio a are analysed. It is found that no matter what kind of jour- 
nal bearings is used, the position of crack ridge zones and then the position of crack unstable zones depend on- 
ly on the speed ratio and the stiffness change ratio for large gravity parameter {W t > 1.0), and depend also on 
the mass ratio for small gravity parameter (IF, <0.05). The position and extent of the crack ridge zones are 
given by Eq.21 for W t > 1.0 and by Eq.22 for IF, <0.05. When 0.1 < 1.0, there is a region, where no 

unstable zones caused by rotor crack exist. 

The larger the speed ratio and the stiffness change ratio are, the wider the crack ridge zone is. Therefore, 
the instability caused by crack in larger speed ratio is more dangerous 

Outside the crack ridge zones, the stability and stability degree of the cracked system arc almost the same 
with those of the uncrackcd system, that is, they have almost no relation with the crack and depend mainly on 
system's parameters; while within the crack ridge zones, they depend mainly on crack stiffness change ratio, 
but also on system's parameters. In some cases, the system may still be stable even the crack is very large. 
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Fijj.l Schematic diagram of a cracked flexible rotor supported on journal bearings 
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(L/D=0. 5) 


FIg.2 Bearing geometry 
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(c) tilting pad, 5 pads 




CilTITT HMITM* ** 

. Cylindrical , 

4 lobti 

(d) unstable limits 


PIr. 3 Influence of gravity parameter on the diagram of stability degree, without 
crack case (AK^ = 0.0, AK^ = 0.0, D« = 0.01 , S 0 « 2.0, a - 0.2) 
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(a) W, =0.001 



SPtl > UTio qm 

(C) Wg = 5.0 


(b) W g =0.1 



(d) unstable limits 


Fig. 4 4-Lobes bearing, influence of gravity parameter and fixed Sommerfcld 
number on the diagram of stability degree, without crack case (AK < = 0.0, 
AK, = 0.0, D, = 0.01 , a = 0.2) 
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Fig.5 Position of crack ridge zones, influence of mass ratio for large gravity 
parameter (4-!obes bearing, AK ( = 0.7, AK,=0.0, D„ = 0.01, S o = 0.1, 
W f = 5.0) 



(a) diagraa of stability degree 
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Fir. 6 Position of crack ridge zones, influence of mass ratio for small gravity 
parameter (Tilting pad bearing, AK { = 0.5, AK„ = 0.0, D e = 0.01, S o = 0.1, 
W r = 0.001) 


320 




***** tf 

(c) Tilting pad, 5 pad* 
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Fig.8 Influence of gravity parameter on the diagram of stability degree 
(AK { = 0.7, AK, = 0.0, D, = 0.01, S 0 = 2.0, a = 0.2) 
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(d) unstable limits 


Fig. 10 Tilting pad bearing, influence of stUTncss change ratio 

parameter on the diagram of stability degree (AK,=0.0, D,= 0. , 

S o =2.0, a = 0.2) 
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